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A COMPUTER PROGRAM 10 CALCULATE THE RfSISlIVIlY Of A THIN FILM OEPOSITfD 
ON A CONDUCTIVE SUBSTRATE FROM F OUR- POINT PROBE MEASUREMENT^ 


Lawrence G. Oberle and Gustave C. Frallck 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


SUMMARY 

This paper deals with the use of the four-point probe to measure the 
resistivity of a thin film of conducting material deposited on another layer 
of conducting material. Such measurements occur, for Instance, In silicon 
carbide (SIC) research, where It Is necessary to grow the SIC on a silicon (SI) 
substrate. The presence of the silicon substrate will Introduce errors In the 
measured resistivity of the SIC. 

Starting from basic principles, and expression for the ratio of measured 
voltage difference to Injected current [AV/I] Is developed. This expression 
Involves the probe spacing, relative thicknesses of the layers, and the sub- 
strate resistivity as parameters, as well as the unknown resistivity of the 
deposited layer. The unknown resistivity can be found by Iteratively evaluat- 
ing the theoretical expression. This must be done numerically. A full 
description of the numerical techniques Involved, and the computer programs 
used, Is given. 

Finally a comparison with previously published results Is presented, 
together with a detailed description of how to use the programs to ' ,nd resis- 
tivities, as well as plot curves displaying the change In AV/I as a function 
of the thicknesses of he layers, and their resistivities. 


INTROOUC1 ION 

One of the ways In which a semiconductor material Is characterized Is by 
the measurement of Its resistivity. In the development of silicon carbide for 
use as semiconductor material for high temperature applications. It became 
necessary to measure the resistivity of the thin film while It still was 
attached to the silicon upon which It had been grown epitaxially. A method of 
calculating the "true" resistivity of a deposited layer on a substrate of 
finite and different resistivity was discussed by Brown and Jakeman (ref. 1). 
The theory was presented In a cursory fashion, and no mention was made of a 
technique (l.e., and analysis, or series of programs) to determine the resis- 
tivity, or conductivity. Our Intent Is to remedy this defect with a presenta- 
tion of the theory In depth, and a library of programs which calculate the 
correction factors necessary for evaluation of the equation presented by Brown. 

The situation Is as shown In figure 1. Current of magnitude I Is Injected 
at probe A and withdrawn at probe B, and the voltage difference, AV, between 
probes 1 and 2 Is measured. The subject of this paper Is the determination of 
the unknown conductivity, o-j , or resistivity, , of the thin film In terms of 
the known quantltltes AV/1, , up (the thicknesses of the respective layers), 

and op (the conductivity of the substrate). 


THEORY 




Governing Equations and Boundary Conditions 
It Is assumed that the materials obey Ohm’s Law, l.e., 

7 = of (1) 

where J Is the current density In A/m?, o is the conductivity (n-m)~7 and 
E Is the electric field In V/m. 

The current density J satisfies the continuity equation for the electric 
charge In the form 


V 



( 2 ) 


where <p Is the electric charge density. In C/m 3 . For steady currents, there 
Is no accumulation of charge at any point, so d<p/dt Is zero, and 


F rom equation ( 1 ) 


The ejectrlc field 
7 x T -> 0, 


V • J = 0 (3) 

V • (oE ) = 0 (4) 

Is derivable from a potential function V, since 


E = - VV (5) 

Combining equations (4) and (5), 

7 • (Tfv) = 0 

or If a Is constant, 

V? V - 0 (6) 

Thus the potential V satisfies Laplace's equations Inside the two layers. 

The boundary conditions satisfied by V may be found by analogy with electro- 
statics. For a charge-free region, Gauss's law (ref. 2) Is given by 

7 • (e7) = 0 (7) 


where c Is the permittivity. Comparing equation (7) with equation ( 4 ), It 
Is seen that a In the steady current problem plays the role of c In the 
electrostatic case Equation (7) leads to the boundary condition (ref. 3). 


8V , 

V 1 ' V ‘1 ST * ‘2 


2 



where a/dn Indicates the normal derivative at the boundary. Replacing 
by a, our boundary conditions become 


C 


av av, 

V 1 3 V °1 in 3 °2 9n 


(8) 


tquatlons (6) and (8), together with the field about a point current source 
(from appendix B), are sufficient tc find the function V. 


Solution of tquatlons 

The method of attack Is to find the potential due to a single point cur- 
rent source and then use superposition to find the potentials due to the two 
sources at A and B. 


Ihe z - 0 plane Is at the top of the upper plane In figure 1, and the 
positive z axis extends Into the material. The two planes are bounded at 
z = 0, z = w] , and z w] ♦ S> n c° we are considering only the source at 

A for the time being, the problem has azimuthal symmetry; that Is, the poten 
1 1 a 1 ha no angular dependence about the z-axls. It Is therefore convenient 
to use cylindrical coordinates. 

In cylindrical coordinates, laplace's equation V - 0, assumes the form 


I a 
r ar 





(9) 


where r Is the distance from the z-axls and o Is the angular displacement. 
Equation (9) may be solved by separation of variables; that Is V(r ,0,z) Is 
written as 


V( r ,0, z) - R(r) 0(d) Z(z) 
In which case equation (9) assumes the form 


1 

rR 



_1 d ? e 
r 2 e do 2 


i ii 

7 , 2 

jZ 


( 10 ) 


( 11 ) 


lhe left side of equation (11) Is a function of r and e, and the right 
side Is a function of z only. In order for equation (11) to be true for all 
values of r > 0, -*< e <«,-•< z < •», both sides must equal a constant. 
Hence equation (11) Is written 


Li. [ dR 
rR dr r dr 


dj0 
t 2 0 de 2 




1 

Z dz 2 


k 


2 


Iherefore 


Z satisfies the equation 


3 


which has the solution 


1 - E|<e^ *■ F ^ e ' 


(l?; 


The solution for e Is found In a similar fashion. Moving the k ? to the 
left side and multiplying by r?, equation (11) becomes 


r d_ 
R dr 



♦ k 2 r 2 


1 d 2 e ? 


(13) 


so that e satisfies the equation 


die 

de 2 


n 2 e 


If the potential Is to be single valued In e, n must be an Integer, and 

e(0) - C n cos (ne) + D n sin (no) ( 14) 

Ihe above choice of separation constants leads to solutions which are periodic 
In 6 and can be made to vanish at z - t ®>. 

This leaves only the equation for R(r). Multiplying by R, equation (13) 
can be rewritten as 


d [ dR*l I - . 2 2 2l „ n 

r L r sfJ * L k r n J R ■ 0 

This Is Bessel's differential equation, and R Is written 

R(r) = A n J n (kr) ♦ B n Y n (kr) 

where J n and Y n are the Bessel functions of order n of the first and 
second kind. Both are oscillatory functions which vanish are r = •, but only 
J n Is finite at r - 0. In our problem, V(0,e,z) must be < ®, so Y n Is 
excluded. Further, since there Is no angular dependence, n = 0, and our 
solution Is written In terms of J 0 (kr). 

For problems Involving boundary conditions on cylindrical boundaries at 
finite values of r and no angular dependence, the solution V(r,z) Involves 
a sum over discrete values of k. As an example, If V Is to vanish on a 
cylinder of radius a, the solution Inside the cylinder Is 

CD 

V(r.z) *£ A e' k 1 ,z, J (k r) (17) 

1=1 01 

where akf = xj, the 1 th zero of J 0 (x). In this way each term In 
equation (17) vanishes at r = a. 


(15) 

( 16 ) 


Now, as a * «>, ki - x^/a becomes d continuous variable, dnd the 
solution (eq. (17)) becomes 

03 

V( r , z) =/ f(k)e~ k|zl J (kr) dk (18) 

0 0 

This Is the form of the solution to be used In finding the potential Inside the 
two conducting layers shown In figure 1, which are assumed to extend radially 
to Infinity. It can be verified that equation (18) satisfies Laplace's 
equation . 


Application to the Four Point Probe 

The potential from the point current source at A will be found first; 
the total potential due to the source at A and sink at B may then be con 
structed by algebraic addition. 

According to appendix B, and the previous disc sslon, the potential In 
the top laver due Just to the source at point A Is 


where 


s 9_ 

J~2 2 

yr * i 



e kz J o (kr) dk 


Q -- 


I 


(19) 


(? 0 ) 


To Vfl are then added terms of the form (eq. (18)) due to the currents flowing 
along the discontinuities at z = 0, and z = wi . If this were an electro- 
statics problem, these additional terms would be due to charge distributions 
at the discontinuities. Combining equations (18) and (19), the potential In 
the top region Is 


V = J [f 1 (k)e kz ♦ g, (k)e~ kz 1 J (kr) dk ♦ - - a - - (21) 

°L J V777 

* fh<- k2 * 9 i <k>e kz * Qe k2 l j 0 (kr) dk (22) 

Likewise, In the bottom layer the potential Is 

V 2 [ f 2 ( k )e kZ f g 2 (»c)e kz J J o (kr) dk (23) 
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Once the funttlons f|, fp, gj, and gp are determined, the potential every 
where Inside the two layers will he known, although, since we need only the 
potential at z - 0, only t| and g t are needed. 

At 7 0, the first of the boundary conditions, equation IS), reduces to 


(r,0) Q 

3Z 


(24) 


since In the region z < 0, a - 0. At z = o>-| , equation (8) gives 

V 1 ( r ' = V 2 (r,<o 1 ) 

aVi(r,U|) aV 2 (r,u^) 

°l az = a 2 dz 

The conductivity for z > <jj ♦ <jp Is zero, so again 

d V ( r ,u * u ) 


(25) 

(26) 


(27) 


The four conditions (24) to (?.l) are sufficient to determine the four unknown 
functions f^, g-|, fp, and g 2 , which appear In equations ( i? J ) to (23). 

Applying equation (24) to (21), 


? V J 

az 




(kr) dk - - 


A L. 


(r 2 ♦ z 2 ) 3/2 


and, at z ; 0, 


av^r.0) 

az 


co 

/ k [ f l< k > V k) ] J o 


(kr) dk = 0 


Since equation (28) must be true for all r > 0, 

fl(k) 9i (k) = 0 

Note that equation (22) was not used In setting 3V-|(r,0)/az = 0 since 


(28) 


( 29 ) 


t "" 


e" J o (kr) dk Is not differentiable at z - 0. 


b 



The next two boundary conditions (25) and (26) give, respectively 

-kw, -kw. 


f • 


'I 

g | e ► Qe I J„( kr ) dk 




([• 


kw, - ku) 


? e ' * g ? e I J Q ( kr ) dk 


■] 


and 


N 


kw. - kw. kw, 

^e - g^e - Qe | J rt (kr) dk 


'h 


no 

j k [* 


kw, -kw. 


' a 2 ) K |V - V I J 0 ( kr ) dk 


Again, since these equations are to be true for all r, the Integrands must be 
equal, whence 


ku.| -kw 


. kw -kw. -kw 

1 1 - 1 - - 1 -- - Qe 1 


^e ♦ g-|C - f 2 e - g 2 e 


and 


kw. 


-kw. 


kw. 


kw. 


c»-| f 1 e 0,9^ ■ ♦ a } q 0 e 


' 1*1 


2 2 2*2 


■ «, o ( e - S ) 


Finally, applying equations (23) and (2/), 


k(w. +w 2 ) k(w. +w ? ) 

f ? e - g ? e - 0 


Rewriting equation (32) as 


kw. -2kw -kw. 

f 2 e = e 2 g 2 e 


[vl 


equations (30) and (31) become, after substitution, 


kw -kw. -kw 

V ♦ g^ - g 2 e 


-2kw / 

1 4 e 2 I 


-kw. 


-Qe 


(30) 


(31) 


(32) 


(33) 


and 
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'ey 


ktii| k(jj ku>j ?koj^ kiU| 

o,f,e o^e f Op9 ? *’ 1 - p = o,Qe 


Eliminating gp between equations (33) anti (34), 


-k« / 2ko> \ / -2k«\ 

V o^l ♦ e 2 j o 2 \l e ? j 


' [-,(< • • *•’) - ;( 


1 - e 


Upon rearrangement, this becomes 
kt» f ? 


■ 

?k(j 


■ 

-2kw 

(o 1 ♦ Op) ♦ (<> ( 

Op)0 

g,e 

_ ( o, 

Op) » (o 1 ♦ Op)P 


Finally, making use of equation (?9) (that f] = gj), multiplying through 
k« 

by e , and dividing through by the quantity (oj *• o ? ) , 


f,(k) = ", 


X_± 


-2ko 

- 2kti)„ 

1 t Re 

R + e 


t,(k) 


1 Re 


2ksa» ? /s 


2ksw./% -2ksw„/s I 1 

e I R ♦ e I 


Qf(ks) 


where we have made use of the definition 


R r V._ a ? . P 2 P 1 

* Op f>2 * P| 

The quantities p^ * l/o^, and p^ - l/o^ are the resistivities. 


We are Interested In the potential difference between the two probes at 
points 1 and 2 In figure 1, when a current I Is Injected at point A and 


8 


« 


withdrawn at point B. from equations (21) and (29), the potential on the 
planp 7 - 0 due to the current entering at A Is 


<r> 

V( r ,0) = J ♦ 2Q / f(ks)J (kr) dk 

t n U 


Since the potentials add algebraically, the potential at point 1 due to 
source Q A at point A, and a source Qg at point B, separated from 
point A by the distance Is 


’P, * S 


f(ks)J o (ks) dk * £ ♦ Q e f(ks)J o (?ks) dk (37) 


Likewise, the potential at point 2, due to the same sources, Is 


V P 2 ’ 2i * 2Q A 


t(ks)J Q (2ks) dk ♦ — ♦ 2Q b f(ks)J Q (ks) dk (38) 


The measured quantity Is the voltage difference between points 1 and 2 


AV r v n V 

P 1 p 2 


^a c r 

♦ 2(U a Qg) J f(ks) |yks) J 0 (?ks)J dk 


In our case, the sources are equal and opposite In s Ign (Qa = - Qb - Q). and 


- ^ ♦ 4Q / f (ks) £j Q (ks) J q ( 2ks)J dk 


Making the change of variable (t = ks) and defining the dimensionless 
quantities 


X „ V- . W - -? 

R o * w* - 


equation (39) becomes, after using the definition of Q found In equation (20) 


AV 1 
1 " 2*s<j, 


1 r 4 


{ ® t 


V t) J o< 2t > 


)J dt 


with the symbolism 
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I 


(4?) 


g(x) 


• «e_ ( 

. e- < “ x > . 


- 1 


The difference J 0 ( t ) J 0 (2t) which appears In equation (41) may be 

evaluated using an Integral representation for J 0 (t) (ref. 4). 


f cos I 

t sln(e)l de 

(43) 

0 1 

J 



whence 


Q (2t) - ~ J £cos |t sln(e)J cos [ 2 1 s1n(e)j j 


J o (t) J 


Using the Identity 

cos (a) ros(0) 2 si 
equation (44) becomes 

2 


e (44) 




J (t) - J ( 2t) = 

O O n 


and equation (41) becomes 


/ s1n|jt 5ln(e)j s 1 n £^t s1n(e)j 


de (45) 


av _ _j. . 

I " 2*o^ s 


1 • * / «G 

a / -b 

t s1n(e)j s1n[~ 

t s1n(e)J de dt 

*'0 

•'o 




where g(x) Is defined as In equation (42). The double Integral, 
equation (46), Is the basis for the programs described In appendix C; the 
details of which are found In the next section. 


PROGRAMMING 

Solution for the Double Integral 

The Integral of equation (46) cannot be evaluated analytically. However, 
Gaussian quadrature (ref. 5) provides a fast, accurate method of evaluating 
the double Integral numerically. In an "m" point Gaussian quadrature, the 
Integrand Is evaluated at points determined by the roots of an m^h order 
Legendre polynomial, and summed, with predetermined weighting factors, to 
provide an approximation to the Integral over the range of Integration. For 
the purposes of our program It was decided that a 6 point quadrature would 
yield sufficient accuracy. 
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Before the Integral over t can be evaluated, the Integral over e 
must be performed. Therefore, we rewrite equation (46) as 


AV P J_ 
I ’ 2 irs 


4 



where, using equation (45) 


K,(t) - J Q (t) - J Q (2t) 


( 47 ) 


or 


^(t) 


f 


Si 


n -t sln(e) sin ^t s1n(e)J 


de 


(48) 


Because of the change of variable In equation (41), K^(t) Is a function of t 
only, and Is evaluated at the same points, regardless of the value of s, pi , 
P 2 , o»i , or w 2 - The points at which K | ( t ) Is evaluated are determined by 
the zeroes of the Integrand of equation (47). Two factors Influence the evalu- 
ation of this Integral. First, an Integral with upper limit of «® Is dif- 
ficult to Integrate numerically because the algorithm to be used requires the 
upper limit be specified. Second, Gaussian quadrature yl ’ds more accurate 
results If the Integration Is performed over subintervals determined by the 
zeros of the Integrand. Both factors are taken Into account with the declsl- n 
to Integrate between the zeroes of the Integrand, and to sum these results 
until the change In the total Is negligible. 


The problem of solving equation (46) Is then divided Into three parts: 
finding the zeroes of the Integrand 


K ? (t) 


g[j] «,<t) 


<«> 


evaluating 


and finally, solving 



(50) 


(51) 


for pi , where Pj Is found by Increasing J until |Pj - Pj*)l Is less than 
the resolution of the computer. 


The variable g[t/x], as given by equation (42), 1 r a monotonlcally 
decreasing function of t, and has only the one zero (at t • •). Therefore, 
the zeroes of Kg are Identical to the zeroes of Kj. The first program- 
ming step, then. 1; to determine the zeroes of Ki(t). The number of zeroes 


11 


of K](t) is determined by the relative c as of the other parameters. We 
found that for all cases attempted, 50 zeroes was sufficient to accurately 
evaluate the Integral. Because equation (48) is a function of t only, the 
zeroes of the Integrand are constant, and can be found once and stored in a 
file. With the zeroes of the Integrand known, the points at which K-|(t) Is 
evaluated are known, and the value of K](t) at these points can be found, 
and stored In a file. 


Ihe 

order of programming Is thus: 



PROGRAM 

1 . 

Determine first "J" zeroes of 

K | ( t ) 


ZERO 

2. 

Determine values of Kj(t) at 

evaluai Ion 

points 

BESCAL 

3. 

Plot AV/1 versus pp/pi for 

various s, 

«! , ti> 2 . and f >2 

FPPPLT 

4. 

Determine p] for given AV/I 

, p 2 , <*n , u>2 

i, and s 

POUR 


Ihe listings for these programs, written In F OR 1 RAN - / / (ref. b) for an I8M-PC, 
can be obtained by contacting COSMIC (Ihe Computer Software Management and 
Information Center f I EM No. 14309J). An explanation of the procedure required 
for running these programs Is found 'n appendix C. 


Program Explanation 

The first two programs, 7ER0 and BFSCAL use the subroutine 8ESDIF to 
evaluate equation (48), at a given value of t. BCSDIF uses the 6 point 
Gaussian quadrature over 30 subintervals spanning the range from 0 to ir. 

ZERO uses an interval halving technique to find the first "j" zeroes of the 
difference, J 0 (t) - J 0 ( 2t ) (with 1" arbitrarily recommended to be 50). The 
first zero Is known to occur at t ^ 0. The program stores this Initial value 
In the ASCII file ZER0.DAT, and begins the algorithm to Identify the rest of 

the zeroes. Starting at t - I, with an Interval of 1, ZERO calculates K^(t) 

until the result changes sign At that point the Interval Is , lved, the 
direction of Increment Is reversed, and the program calculates values of K^t) 
until the sign changes again. This process repeats until the desired accuracy 
(a recommended 5 digits) Is achieved. When the zero Is found sufficiently 

accurately, the value of t which produced the zero Is stored, the Interval 

Is reset to 1, and the program repeats until the desired number of zeroes Is 
found. 

Ihe program BESCAl , using the location of these zeroes, then calculates 
the values of K](t) necessary for the other two programs to evaluate the 
Integral of equation (4b). For each Interval between zeroes, BESCAL cal- 
culates and stores the values of Kj(t) In the file TWOBES.DAT. 

With the locations of the zeroes known, and the values of K-j ( t ) known 
at the points of evaluation, the two main programs, FPPPLT and FOUR can be 
executed. FPPPLT calculates an array of numbers suitable for plotting. The 
plots appearing In this report were generated using an off-the-shelf spread 
sheet program (ref. 7). The plots produced are log(AV/I) versus 1og(/>2^Pi) 
for constant W, s and p?, and for various values of x. This program was 
used by the authors primarily as a check of the algorithms used. 
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r 


f 


FOUR, the more useful of the two programs, calculates a to corre 

spond to a given wi , o>p, p?, s, and AV/I. This Is accomplished by multiplying 
both sides of equation (47) by ?« s/p 



(52) 


The left side of equation (52) can be evaluated from experimentally determined 
parameters. The right side of equation (52) Is calculated from an Initial 
guess of p ] , keeping In mind that g(x) Is a function of p 1 , as given by 

equations (42) and (36). Theoretically the right side of equation (52) can be 
re-evaluated repeatedly until the equality Is achieved. In practice the Inter 
mediate functions f' and F are calculated as: 


f ' 



(53) 


and 



(54) 


If | F 1 - F | Is less than the required tolerance, chen the value of 
p 1 used Is displayed as the "correct" resistivity. Otherwise, the guess for 

p 1 Is modified accord’ng to the sign of the quantity F' - F, and If neces 

sary, the Interval Is halved, until the error Is less than the required toler- 
ance. Because most of the "number crunching" needs be performed once for all 
possible values of the variables, the calculation speed of the algorithm Is 
very high. 


RESULTS 

Ihe curves of figure 2 were produced using data generated by the program 
FPPPL1. The parameters were chosen to correspond to the parameters chosen by 
Brown (ref. 1), In order to compare figure 2 with Brown's results. The shapes 
of the curves are Identical to those produced by Brown. The difference In the 
plots Is that the graphs of this paper are of AV/I versus pp/pj , and In 
Brown's report, the plots are of the function F (Identical tc F In this 
paper) versus o^/op. The abscissa Is the same In both cases, but the 
ordinate In Brown's paper does not correspond to a measurable quantity. Since 
F' Is a function of the probe spacing, and the resistivity of the substrate, 
these two parameters were set to one In figures 2(a) to (c). With this choice 
of parameters. 


F 


AV ,, , 
( 2 ») 


(55) 



» 

'» 
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laklng this Into account, the plots of this paper are Identical with Brown's 
plots . 

figure 3 shows the variation of AV/1 as a function of the ratio of the 
resistivities of the layers as well as a function of the thickness of the 
deposited layer. When the resistivities are approximately equal (pp/pi 3 1 ) . 
AV/I Is a function of the thickness of the deposited layer only In the sense 
that AV/I Is dependent on the thickness of the sum of the two layers 
(w] * up) • Intuitively, this can be seen to be correct, for If the resis- 
tivities are equal, the division between the materials disappears, and the 
measurement of AV/I Is that of a single layer. Since the deposited layer 
varies from less than 2 percent to less than 10 percent of the total thickness, 
the variation In AV/I is relatively small. As the ratio p£/p^ Increases, 
the thickness of the deposited layer comes more Into play. In general. It can 
be seen that the thicker the deposited layer, the lower the measured value of 
AV/I. In addition, the slope of the curve steepens as the ratio p;>/pi 
I ncreases. This Implies that the algorithm Is more accurate for values of 
PI <<: P2- 

In order to verify the algorithm In the program FOUR, the two examples 
from Brown's paper were chosen to be used In our calculation. Table I 
reproduces the applicable numbers from Brown's table II, augmented by the 
results obtained using the program FOUR. 

Ihe results of our program agree closely with the results published by 
Brown, and with the known values for the substances In question. The slight 
discrepancies are caused by the Increased resolution of our machine, and by 
assumed roundoff errors associated with the numbers published by Brown. With 
out knowing the algorithm used by Brown, it Is difficult to make Judgments 
about his techniques. However, It must be noted that advances In computer 
technology, In the 20 years since Brown reported his results, have made the 
evaluation of these equations much easier, and more accurate than before. In 
addition to this Is the fact that the program FOUR solves the equation for p, 
and Brown's results are given In terms of 1/p = a. 


CONCLUSIONS 

As can be seen from the plots, and from the numbers generated to produce 
table 1, the algorithm devised solves the problem as well as the unknown algo 
rlthm designed by Brown. Along with a detailed explanation of the software 
generated In this effort, we tried to explain the process of arriving at the 
Integral presented by Brown, starting from the governing equation, and boundary 
conditions. It Is hoped that the details In the solution are sufficient to 
enable the reader to use these results In the solution of the class of problems 
where resistivity measurements are needed for a substance deposited on a sub 
strate of higher resistivity than the deposited layer. 


APPENDIX A 
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An 
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Bn 

^n 


On 

r 

tk 

F 

F 1 

fk 

f 

h 

<2 

9 

91 

92 
I 

7 

Jn 

Jo 
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NOMENCLATURE 
probe at which current Is Injected 

coefficient of 1st Independent solution to Bessel's equation, V 
probe at which current Is withdrawn 

coefficient of 2nd Independent solution to Bessel's equation, V 

coefficient of 1st Independent solution for 0(G) 

constant to be determined In solving for potential due to an 
Isolated point source 

coefficient of 2nd Independent solution for 0(G) 
electric field, V/m 

coefficient of 1st Independent solution for Z(z) 

Intermediate function used In computing resistivity 

Intermediate function calculated from measured parameters 

coefficient of 2nd Independent solution for Z(z) 

function to be determined In solving the potential problem In the 
continuous case 

function solved In the deposited layer due to current Input at 
point A 

function solved In the substrate due to current Input at point A 

function f redefined In terms of x rather than ks 

function (akin to f^) In the deposited layer due to current 

withdrawn at B 

function (akin to f p) In the substrate due to current withdrawn 
at B 

current Injected at probe A, A 
current density, A/m ? 

Bessel function of the first kind of order n 

Bessel function of the first kind of order 0 

number of zeroes of K](t) to be found 
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A 


K 1 

k 2 

k 

m 

n 

^m 

Q 


Qa 

Qb 

R 

R(r) 

r 


V(r,e,z) 

Va 

V 


'1 


'1 


V? 

w 

Y n 

2(0 

7 

10 


difference J 0 (t) - J 0 (?t) 
g(t/x) k t ( t ) 

separation constant, m- 1 

order of legendre polynomial used In Gaussian quadrature 

separation ronstant, nr ^ 

value of Im ^2 dt using "m" zeroes 


Psuedo charge (analogous to charge In electrostatics) 
Q = I/( 2iro-| ) , V/m 


Psuedo charge at point A, V/m 
Psuedo charge at point B, V/m 


resistivity variable ( P? Pj)/(p 2 * 


potential function In the r direction, V 

radial coordinate 

unit vector In r direction 

spacing between any two adjacent probes In four point apparatus, m 

elemental area on surface of sphere 

potential function, V 

potential function at probe A, V 

potential function at point 1, V 

potential function at point 2, V 

potential function In deposited layer, V 

potential function In substrate, V 

ratio of substrate thickness to thickness of deposited layer 

Bessel function of the second kind of order n 

potential function In the i direction 

axial coordinate 

absolute value of i 
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'N . 1 




AV 


‘I 

c 2 

0 ( 0 ) 

e 

x 

p 

p> 


<p 


w 2 

a 

T 

v 2 


voltage cl 1 f t front. e between probes 1 and ?, V 
permittivity, f/m 
permittivity of first layer, 1 /m 
permittivity of second layer, f/m 
potential function In the ©direction 
angular coordinate 

twice the ratio of probe spacing to deposited layer thickness 
resistivity, ft-m 

resistivity In the deposited layer, Q-m 
resistivity In the substrate, Qm 

conductivity, (tim) •, sr W 1 
conductivity of the depc ited layer, 

conductivity of the substrate, 

electric charge density, C/m^ 
thickness of the deposited layer, m 
thickness of the substrate, m 
partial derivative operator 
vector differential operator 

Laplaclan operator 


17 


k 


APPENDIX B 


Ihe Electric Potential due to an Isolated Point Current Source 


Consider a point current source of output I located at an origin within a 
medium of uniform conductivity, a. By symmetry, the current density vector, 

J must vary as l/r?, where r Is the distance from the source, since the 
total current passing through any sphere surrounding the source Is I, and the 
area of the sphere Is proportional to r ? . That Is 


• ds = I 
sphere 


where ds Is the element of area on the surface of the ,phere, 


(B.l) 


- - 2 

ds = r r sin e de dtp 


( B . 2) 


Since J = ot, the electric field E Is also directed along r and Is 

2 

proportional to 1/r . Thus 


t - 


c 


2 


r 


The unknown constant, c, In equation (B.3) can be determined from 
equations (B.l) and (B.2): 



sin e do dtp 


- oc 



sin 0 do dtp = 4 *oC 


(B.3) 


(B.4) 


or 


c 


_ !_ 

4wo 


Then, from equation (B.3), the electric field becomes: 


( B. 5) 


E 


. 1 

4*o 


* 

r 



( 8 . 6 ) 


Ihls Is the electric field about an Isolated current sorrce of strength I. 
Since the potential, V, Is related to T by T * - ?V, the potential 
(referenced to zero at Infinity) Is: 


f 



. A 1 / J >Vlk. * l - 



1 he electric field, as given by equation (B.6), and the corresponding 
potential, from equation (B.7), apply to the case of an Isolated point source 
from which the current density Is Independent of direction. In our case, 
however, the current (In cross section) Is confined to the upper half plane 
(o = 0 for z < 0) . Then the limits on e In equation (B.4) are 0 and ir/2, 
and equation (B.5) becomes: 


and the potential In the top layer due just to the current Injected at point 
A Is: 


In cylindrical coordinates, where r Is the distance from the z-axls, 
equation (B.9) becomes: 


V 






I 

2ir o 


tx 

L 


e~ kZ J Q (kr) dk 


(B.10) 
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APPENDIX C 


Operation of Programs and Sample Outputs 

It Is assumed that the programs are all compiled using the Microsoft 
FORTRAN- 77 compiler on an IBM PC or compatible computer. The executable 
programs are then generated by using the Microsoft Linker to link the object 
code In the following manner: 



PROGRAM 

1 INKED OBJECT MODUL 

1 . 

ZERO 

ZERO, BESOIF 

2. 

BESCAL 

BESCAL, bESDIF 

3. 

FPPPLT 

FPPPL1, QUAD, XI 

4. 

FOUR 

FOUR, QUAD, XI 


lo run these programs on an IBM PC simply enter the name of the program 
to be run. ZERO and BESCAL need only be run once to generate the reqtlred 
constants for the other two programs. The following pages show example outputs 
of the program as run on an IBM PC. 

The first program to be run Is the program ZERO. The prompts are: 

ENTER REQUIRED NUMBER OF 0IG1IS [DEF = 50] : 

(Your response: an Integer greater than 0.) 

ENTER REQUIRED DIGITS [DEF = 5]: (accuracy) 

(Your response: an Integer between 1 and 6, Inclusive.) 

The program proceeds to display the number of zeroes requested, as well 
as storing them In the file "ZERO. DAT". 

lhe next program Is BESCAL. There are no prompts, but the program d 1 s 
plays the Interval for which It Is presently computing the necessary values 
of K^(t). These values are stored In the file "TW0BES.DAT". 

lhe third program, FOUR Is the program which will be most useful to the 
researcher. The prompts which appear In this program are: 

INITIALIZING QUADRATURE VALUES (reading the files "ZER0.DAT" and 
"TW0BES.DAT" . ) 

Enter PROBE SPACING <0.159>, cm: 

Enter MEASURED VOLTAGE/CURRENI <36.7/>, 0: 

Enter THICKNESS of SI <381. 00>. urn: 

Enter THICKNESS of SIC <6.00>. v m: 

Enter RESISTIVITY of SIC <0.22 E ♦ 02>, Q-cm: 

Enter INITIAL GUESS for RESISTIVITY OF SIC <0.1 E ♦ 01>, B-cm: 

Enter INITIAL DELTA RESIST I VITY <0.1 E ♦ 01>, fl-cm: 

Enter MAXIMUM No. of ITERATIONS <100>: 

CHANGES [Y/N] <N>: (If changes are requested the rroQram loops through 

these prompts again.) 
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lhe routines then 

produces 

a table 

similar to 

the one bel 

1 

Fc 

Fe 

ERR 

RH01 

1 1 

.30009 

1.6697 

0.22138 

0. 1000E+00 

2 3 

.69198 

1.6697 

1.21111 

.600E+00 


11 

1 .59509 

1 .6697 

.04470 

. 1 31 3E+G0 

12 

1.66295 

1 .669 7 

.00406 

. 1 391 E+00 


The calculated SIC RESIST I VI 1 Y = 0.13906E*00 G-cm. 

Following the solution appears the prompt: 

ANOTHER MEASUREMENT ? [V/N] <Y>: 

(If the response Is "Y" the program returns to the original prompts; If 
the response Is "N" the program ends.) 

The fourth program In the package, FPPPLT, Is similar In operation to the 
program FOUR. This routine produces a file which can be used by a package such 
as "SYMPHONY" to produced plots of the relationship described In this paper, 
lhe prompts are: 

INITIALIZING QUADRATURE VALUES (reading the files "ZER0.DAT" and 
"TWOBES.DAT".) 

Enter No. of INCREMENTS <20>: (number of points to plot) 

Enter PROBE SPACING <1.00>, cm: 

Enter THICKNESS of SUBSTRATE <300. 00>, urn: 

Enter THICKNESS of DEPOSITED LAYER <20.00>, pm: 

Enter RESISTIVITY of SUBSTRATE <0.10 E f01>, G-cm: 

Enter BEGINNING 1og(RH01/RH02) <0.00>: 

Enter ENDING Tog(RH01/RH02) < 3.00>: 

CHANGES [Y/N] <N>: (If changes are requested the program loops through 

these prompts again.) 

At this point the routine requires a file name under which to store the 
data. The prompt for this Is: 

ENTER FOR UNIT 2 THE NAME OF THE OUTPUT FILE 
WHICH WILL ACCEPT THE DATA TO BE PLOTTED. 

File name missing or blank - please enter name. 

UNI I 2? (Your response Is to be a legal file name, such as "TEST. OAT".) 
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TABLE 1. - COMPARISON OF THEORY WI1H EXPERIMENT 


FOR A BRASS SHIM ON A MERCURY SUBSTRATE 


S, 

cm 

X 

W 

°2» 

k/(Q- cm) 

AV/I, 

pV/A 

Known 

value 

[o^ - k/(Q-cm)] 
Derived 
BROWN 

Derived 

FOUR 

0.127 

0.063 

4.72 

2.27 

23.2 

19.0 

10.4 

10.4 

63.0 

82.0 

143.0 

141.0 

149.0 

137.0 

150.0 

138.0 


Note: k/(ft-cm); denotes thousands of (Q-cm)" 1 , l.e., 

1 k/(Q-cm) - 103 (Q-cmH 


Also: For a * 10 k/(Q-cm) 

p » l/o * 10-*, (Q-cm) 
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ORIGINAL PAGE !3 

OF POOR QUALITY 




Figure 1. - Theoretical setup for four-point probe measurement of deposited layer on 
substrate of finite conductivity. 
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